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ABSTRACT 

Detailed heat transfer measurements and predictions '/ 

are given for a turbine rotor with 13()° of turning and r 

an axial chord of 12.7 cm. Data were obtained for inlet R ( in 

Reynolds numbers of 0.5 and 1 .0 x 10*5 for isentropic exit ,s 

Mach numbers of 1.0 and 1.3. and for inlet turbulence M 

intensities of 0.25% and 7.0%. Measurements were made f 

in a linear cascade having a highly three-dimensional flow ^ 

field resulting from thick inlet boundary layers. The pur- ^ n 

pose of the work is to provide benchmark quality data ^ 

for three-dimensional CFD code and model verification. V 

Data were obtained by a steady-state technique using a. r 

heated, isothermal blade. Pleat, fluxes were determined 7 

from a calibrated resistance layer in conjunction with a ^ 

surface temperature measured by calibrated liquid ervs- ( 

t.als. The results show the effects of strong secondary ^ 

vortical flows, laminar-to-turbulent transition, shock im- 
pingement. and increased inlet turbulence on the surface P 

lieat transfer. P* 

P 

(T 

LIST OF SYMBOLS Subscripts 


Cp 

specific heat [J/kg-K] 

aw 

C, 

blade axial chord [cm] 

Cu 

d 

rotor leading edge diameter [cm] 

Ic 

k 

thermal conductivity [W/ni-K] 

e.r 

k+ 

equivalent roughness height 

2 7} 

L.r 

longitudinal integral turbulence length scale [cm] 

IS 

M 

Mach number 

Superscripts 

P 

pressure [Fa] 

/ 

Pr 

Prandtl number 



heat flux [YY/ni"] 
recovery factor, r — Pr 1 /’* 

Reynolds number. R< m V/p 

blade surface coordinate [cm] 

Stanton number 

thickness [cm] 

temperature [K] 

turbulence intensity 

total velocity [m/s] 

equivalent normal distance 

spanwise (radial) coordinate 

specific heat ratio, ■*)• — I A 

99% boundary layer thickness [cm] 

emissivity 

nondimensional surface temperature 
= (T lc -T aw )/(Tcu-T {c ) 
dynamic viscosity [kg/s-m] 
turbulent eddy viscosity [kg/s-m] 
density [kg/nv 1 ] 

Stefan-Boltzmann constant 

adiabatic wall temperature 
copper substrate 
liquid crystal 
exit freestream value 
inlet freestream value 
isentro]>ic value 

total conditions 
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INTRODUCTION 

Highly detailed heat transfer data at conditions ap- 
proximating those in actual engines are needed to ver- 
ify computational fluid dynamic (CFD) predictive analy- 
ses. Accurate local heat transfer predictions are required 
to improve blade service life and to reduce cooling re- 
quirements. Accurate heat transfer measurements for 
CFD verification require the good spatial resolution of 
a large-scale facility. The large scale, combined with high 
Reynolds and Mach numbers, in turn requires high mass 
flow rates per passage. A linear cascade, with fewer blades 
t han a full annular cascade, provides better spatial resolu- 
tion for the same overall How rate. Previous studies have 
shown that rotor geometries in linear cascades provide 
good midspan data as compared to their rotating equiva- 
lents. See, for example, Bauglin et al. (1995) or Guenette 
et al. (1989). Also, Graziani et, al. (1980) studied the ef- 
fects of thic k and t hin endwall boundary layers on rotor 
heat transfer in a low speed linear cascade. Blair (1994) 
measured heat transfer in a low speed facility on a rotat- 
ing blade that had the same midspan section as the blade 
of Graziani et al. (1980). Similar three-dimensional heat 
t ransfer patterns were observed at and below midspan of 
those two tests. Linear cascades can and have also been 
used to st udy tip clearance heat t ransfer, but the current 
study emphasizes blade data away from the tip region. 

Although many turbine applications use film cooling 
for blade temperature control, data without film cooling 
are needed to validate predictions of heat transfer coeffi- 
cients. With film cooling present , differences between the 
analysis and data could be due to differences in either 
lirat- transfer coefficients or in the film effectiveness dis- 
tribution. Transition influences blade heat transfer, while 
endwall flows are more likely to be turbulent. Transition 
is not an issue in film cooled blades, where the film cool- 
ing trips the flow to turbulent. However, it is important 
to accurately predict transition, both for heat transfer on 
a non-film cooled blade, and for the aerodynamics of the 
low pressure turbine where heat transfer data can be used 
to verify transition predictions. Secondary flow effects on 
heat transfer differ between the blade and endwall. 

Aerodynamic' and heat transfer data are needed at 
the same flow conditions to thoroughly understand the 
flow physics. The objective of the current study is to 
provide a detailed, high Mach number rotor blade heat 
transfer data set to CFD code developers and users so 
that they can better manage the uncertainty associated 
with heat transfer predictions. The current study, along 
with the aerodynamic study (Giel et ah, 199<>a), and the 
endwall heat transfer study (Giel et al., 1998), provide a 
complete and detailed data set for CFD code and model 
validation. A baseline 3-D CFD analysis is included to 
illustrate where improved modeling is needed. 


Along with the studies mentioned above, de- 
tailed measurements were reported by Goldstein and 
Spores (1988) and Goldstein et a). (1995) for rotors in 
large scale low speed linear cascades. These studies 
used a naphthalene sublimation technique while Graziani 
et al. (1980) used heated blade and endwall surfaces. 
Dunn et al. (1994) measured rotor heat transfer at dis- 
crete chordwise and spanwise locations in a shock tube fa- 
cility at. engine-typical gas- to- wall temperature ratios and 
Mach numbers. These measurements were on engine-size 
hardware using heat flux gauges. Consequently, the res- 
olution was less than for the large scale facility measure- 
ments. Martinez-Botas et al. (1994) used a liquid crystal 
technique to measure heat transfer for an uncooled stator 
in a blowdown annular cascade at transonic conditions. 

The NASA Glenn Research Center Transonic Turbine 
Blade Cascade provides heat transfer and aerodynamic 
data to verify CFD analyses. Data were obt ained for eight 
different flow conditions, and illustrate the effect s of vary- 
ing Reynolds number, exit Mach number, and inlet turbu- 
lence. Measurements were made at inlet Reynolds num- 
bers of 0.5 and 1.0 x 10'\ and exit pressures correspond- 
ing to isentropic Mach numbers of 1.0 and 1 .3. Tests were 
conducted at high and low inlet turbulence levels. Twelve 
blades ensured good periodicity. The axial chord was 
12.7 cm. The Reynolds and Mach numbers are relevant 
to high pressure turbines where heat transfer is critical. 
Results herein show the effects of shock/boundary layer 
interactions and of longitudinal vortices on blade heat 
transfer, which become more significant as Mach num- 
bers and flow turning increase. 

The rotor, of constant cross section, is highly loaded 
with 136 degrees of turning, and an inlet flow angle of 63.6 
degrees. The cascade geometry is two-dimensional, but 
the flows are highly three-dimensional due to thick bound- 
ary layers developed on the long cascade inlet. Each 
boundary layer was one-third of the half-span. 1 he full 
span-to-axial chord ratio was 1.2. The thick boundary 
layers and high flow turning resulted in spanwise flow vari- 
ations as large as those seen in rotat ing turbines. See, for 
example, Joslyn and Dring (1992) or Thulin et al. ( 1982). 
These variations result from secondary flows, and accu- 
rately predicting them is a significant test of a three- 
dimensional analysis. Tests were done with and without 
a blown, square bar grid upstream of the blade row. 

Using liquid crystals to measure heat transfer gave 
good spatial resolution. Several versions of liquid crystal 
measurement techniques have been used for gas turbine 
related measurements. A transient technique was used 
by Martinez-Botas et al. (1994), and a steady-state tech- 
nique was used by Hippensteele et al. (1985). The rough- 
ness of the liquid crystals is documented because of its 
heat transfer importance. 
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Fig. 1 Transonic Turbine Blade Cascade 

Rotor shape was influenced by both aerodynamic 
and heat transfer considerations. Predicted rotor profile 
loss was consistent with current design practice. Arts et 
al. (1997), Dunn et ah (1994). and Blair (1994) among 
others, showed midspan rotor blade leading edge heat 
transfer rates nearly twice the highest values seen else- 
where on the blade. Civinskas et al. (1990) showed heat 
transfer in this region was reduced using a large leading 
edge diameter, and aerodynamic losses were not signifi- 
cantly increased. The rotor shape was chosen to achieve 
acceptable aerodynamics with a lower than typical ra- 
tio of peak-to-average heat transfer. This lower peak- 
to- aver age heat transfer had a secondary benefit in that 
measuring heat transfer over a narrower range inherently 
reduces experimental uncertainty. 

DESCRIPTION OF FACILITY 
Cascade Description An initial facility description 
was given by Verhoff et al. (1992). The cascade inlet 
section was analyzed and redesigned ((del et al. 1990b) 
to improve inlet flow pitchwise uniformity. Figure 1 shows 
an overall view of the redesigned facility. Passage I is at 
the upper left and passage 11 is at the lower right. 

High pressure, ambient temperature air was throttled 
to an inlet total pressure near one atmosphere for the high 
Reynolds number tests. Low Reynolds number tests had 
an inlet tot al pressure near one-half atmosphere. The air 
discharged into an exhaust header maintained at a nom- 
inal pressure of 15.9 kPa (2.3 psia). A valve between the 
test section and exhaust header was used to maintain the 
desired exit Mach number. Figure 1 shows the test sec- 
tion mounted on a large disk, which can rotate to vary 
inlet incidence angle. Results herein were obtained at the 
design inlet flow angle of 63.6 degrees. Upstream inlet 
boards were used, but to minimize downstream shock re- 
flections no exit tailboards were used. Thick boundary 
layers developed in the long inlet section. Aerodynamic 
probe and blade loading data verified the existence of 
strongly three-dimensional passage flow. 
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Fig. 2 Heat Transfer Measurement Blades 


Table 1 Blade and cascade parameters 
and dimensions 


(Geometric parameter 

Value 

axial chord 

12.70 cm (5.000 inches) 

pitch 

1 3.00 cm (5. 1 1 9 inches) 

span 

15.24 cm (0.000 inches) 

true chord 

18.42 cm (7.250 inches) 

stagger angle 

4 1 .54° 

throat dimension 

3.358 cm (1.393 inches) 

throat area: l passage 

53.94 cirr 

leading edge diameter 

2.057 cm (1.040 inches) 

trailing edge diameter 

0.518 cm (0.204 inches) 

turbulence grid 

2.54 cm square bar 

Flow parameter 

Value 

Inlet Reynolds No. 

0.5 x 1() G or 1.0 x 10° 

Exit Reynolds No. 

0.9 x 10 ,; or 1.8 x ){)" 

Inlet Mach No., Mj$ 

0.38 

Exit Mach No., Mis 

1.32 or 0.98 

Inlet S - no grid 

3.2 cm ( 1 .2 inch) 

with grid 

2.0 cm (0.8 inch) 

Inlet flow angle 

03.0° 

Design flow turning 

130° 


Blade and cascade details are given in Table 1. The 
two blades forming passage 5 were instrumented; one for 
the suction surface and one for the pressure surface with 
some leading edge overlap. Both endwalls were 0.35 cm 
(2.50 in.) thick clear acrylic for optical access. 
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Measurement Blade Description High strength 
metal blades were fabricated from high-conductivity, 
oxygen- free copper to withstand large pressure differen- 
tials due to high exit Mach numbers. They were under- 
cut by 0.(56 mm (0.02b in.) to accommodate a. composite 
low thermal conductivity layer across which heat trans- 
fer rates were measured. Figure 2 shows the layer to be 
0.15 mm (0.000 in.) of double-faced adhesive film and 
0.51 mm (0.020 in.) of neoprene rubber. Liquid crystals 
were sprayed on the outer surface. The adhesive film and 
the neoprene rubber were rolled carefully onto the cop- 
per surface. The following transient tests were performed 
to verify uniform bonding of both: A heat gun was used 
to quickly heat, the surface until the liquid crystal yellow 
line was visible. The surface was then allowed to cool. 
Any regions with air bubbles trapped in the composite 
layer responded more slowly to the transient heating and 
cooling and thus appeared as obvious distortions in the 
yellow line. The layer in this region was then punctured 
with a fine hypodermic needle to eliminate the bubble. In 
another test , the blade heaters were turned on abruptly 
while the test section was under a vacuum. A spanwise 
liquid crystal yellow line was observed to move from the 
thinner sections of the blade slowly towards the thicker 
sections. The two-dimensionality of this line further ver- 
ified uniform composite layer bonding. 

Cylindrical electrical cartridge heaters extended the 
entire span of the blades just under the surface. They 
were positioned to maintain copper temperature unifor- 
mity within 0.1 5°(' 4 (0.25°F). Sheathed Type E 0.51 mm 
(0.020 in.) thermocouples, ill copper surface grooves, 
were used by feedback control circuits and a D.C. power 
supply to maintain the copper substrate at a uniform tem- 
perature. For most flow conditions the midspan copper 
thermocouple temperatures were maintained to within 
±0.6°C (±1°F). An additional row of thermocouples at 
t he same surface distances as the midspan row was lo- 
cated (i.4 mm (0.25 in.) from the end wall nearest the 
cameras. The four thermocouples surrounding any point 
on the blade were used in a bilinear (s, c) interpolation 
procedure. 'This interpolation gave the local Tc u and thus 
the local heat flux at any surface location, eliminating the 
need to correct for end wall heat loss. 

The composite layer conductance, k/t , was calibrated 
in a separate test. A small rectangular test plate was 
made of the same copper material and was instrumented 
with the same thermocouples, adhesive film, neoprene 
rubber, and liquid crystals that were used in the blade 
tests. A thin-fihn electrical heater was bonded to the 
bottom surface and the current and voltage were care- 
fully monitored to determine the heat flux through the 
plate. Thermal insulation, thermocouples, and guard 
heaters were placed on the bottom surface and around 


all four edges. Compressed air was blown on the sur- 
face at various flow rates to in order to vary the external 
heat transfer coefficient . At each blowing rate, the heater 
power was adjusted until the liquid crystal yellow line 
was visible. The conductance was then calculated from 
the heat flux and the temperature difference between the 
copper substrate and the liquid crystals. The composite 
layer conductance. A*//, was calibrated as 4(58 W/m J *K, 
±5% (82.5 BTIJ/hr- ft- R) and was not found to vary 
significantly over the range of heat fluxes tested. The 
heat transfer measurement technique was first tested on 
a circular cylinder in cross-flow, and results agreed with 
well-established data. 

A profilometer measured several 5-mm traces after 
the liquid crystals were sprayed on. A (j(Jx photomicro- 
graph showed that the profilometer stylus did not alter 
the crystals. The arithmetic mean roughness was (5.5 /mi 
and the r.m.s. roughness was 7.8 //m. The average spa- 
tial frequency of roughness was 10.4 mm' so the aver- 
age peak-to-peak distance was 0.097 mm. The equivalent 
roughness height was estimated to be 7 /an. The maxi- 
mum normalized roughness height, A + , was estimated to 
be 3.5, so from a surface roughness standpoint the blade 
surface was hydraulically smooth. More details of the sur- 
face characteristics including digitized traces and surface 
FFT's are available from the authors. 

Turbulence Grid Description A 74% open area 
square bar turbulence grid. (5.91 upstream of the mea- 
surement passage, positioned normal to the inlet flow, 
was used for some tests (see Fig. 1). Square 25 mm 
(1 in.) tubes, one extending between the inlet boards 
at midspan, interconnected to three spanwise tubes lo- 
cated on 150 mm centers formed the grid. Balanced air 
at 965 kPa (125 psig) entered the entire grid through 
both ends of the spanwise tubes, and exited upstream 
through 75 3.2 mm (1/8 inch) diameter holes. The holes 
were spaced uniformly on 12.7 mm (1/2 inch) centers 
with 11 on each spanwise tube and 42 on the tube ex- 
tending between the inlet boards. Each row of holes was 
centered on its respective square tube. The flow from 
each hole directly opposed the mainstream flow. Total 
mass flow from the grid was 0.(55 kg/s (1.4 lb ni /s) for all 
cases, corresponding to 5% of the cascade mass flow at 
R( i n = 1.0 x 10 6 . This grid was used in a previous end- 
wall heat transfer study (see Giel et, ah, 1998), but with 
no grid air. Boyle et al. (1998a) showed that upstream 
blowing produced a more uniform mean and fluctuating 
flow field compared to an unblown grid. 

MEASUREMENT TECHNIQUES 
Inlet Flow Measurements Aerodynamic probe mea- 
surements were made on a plane one axial chord up- 
stream of the blade leading edge plane. In the pitch- 
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wise direction, the measurements covered the primary 
heat transfer measurement passage, passage 5, and half 
of the adjacent passages 4 and (5. Span wise, the mea- 
surements started near the endwall, c /span ~ 0.042, and 
extended to midspan, z japan = 0.50. Time mean flow 
measurements were made with a calibrated 5-hole pitch- 
vaw probe. Details of the measurement techniques are 
given by (riel et al. (1990a). Turbulence intensity and 
integral length scale were measured with a single 5 //in 
hot wire. The wire output voltage was calibrated against 
a density-velocity product in a free air jet; this compen- 
sated for the fact that pressure in the rig was lower than 
atmospheric. Data was recorded at a sampling frequency 
of 70.2 kHz with an antialiasing filter cutoff frequenc y of 
25.4 kHz. 1048570 data points were recorded for each 
length scale. The data was analyzed in 32768 blocks with 
50% overlap. This resulted in a lowest distinguishable 
frequency of 2.3 Hz. The resulting 03 autocorrelations 
were then averaged and the result fit with an exponential 
curve of the form r(i) — exp (—Ct). Length scale was 
then computed as L r = ( ' aV g/(. 1 as discussed in Van Fos- 
se n et al. (1995). 

Blade Static Pressure Measurements Surface 

static pressures were measured using a separate blade set. 
The t wo blades forming passage 5 were instrument ed wit h 
149 stat ic pressure taps on 9 spanwise planes extending 
over the entire blade span. Details of the blade static 
pressure measurements are given by (4iel et al. (1990a). 
Heat Transfer Measurements An array of dots, vis- 
ible through the liquid crystals, were first, painted on the 
outer surface. Next, two micro-encapsulated chiral ne- 
matic' liquid crystals with different vellow-band tempera- 
tures were mixed and sprayed onto the surfaces, and onto 
a separate instrumented calibration plate to determine 
their temperature characteristics. Moffat (1990) showed 
that micro-encapsulated crystals are less sensitive to illu- 
minating and viewing angles than non-encapsulated crys- 
tals. No differences in isotherm location were detected 
between overlapping photographs. 

Two liquid crystals, each with a full color bandwidth 
of ±1°C, were used due to surface temperature limits. 
Low surface temperatures give high uncertainty, and high 
surface temperatures cause the neoprene to separate from 
the surface. The higher temperature crystal measured low 
Stanton number regions and the lower temperature crys- 
tal measured high Stanton number regions. For an inlet 
air temperature of 18°C (65° F), crystal temperatures of 
38.1°C (100.6°F), and 48.1°C (118.5°F) were used. In re- 
gions where data were available from both crystals, agree- 
ment between the crystals verified that the correct defini- 
tion of Stanton number was used and that the adiabatic 
wall temperature was reasonably accurate. 
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Fig. 3 Sample digitized liquid crystal data distribution 
(Re jn = 0.5x10 6 ; M ex =1.32; no grid) 

When steady state was achieved, four 35 mm color 
slide cameras photographed the blade surface. Each cam- 
era viewed a subset of the blade surface, with some over- 
lap between camera views. High speed photographic 
strobes were used to prevent radiative heating of the 
liquid crystals. The photographic slides were projected 
onto paper marked with the same dot, pattern as on the 
blades. Slide image dot patterns were aligned with the pa- 
per dot patterns. Drawings of the isotherms were made 
for each camera and crystal temperature. These lines 
were digitized, with between 2000 to 5000 digitized points 
for each flow condition. Each point was mapped onto 
the unwrapped blade coordinates, and the mapping cor- 
rected for distortions due to blade curvature. Typically, 
20 heater sett ings covered the full range of surface Stan- 
ton numbers. Figure 3 shows a sample digitized data dis- 
tribution. (4 ray areas are regions with no available data 
because of shadowing from other surfaces. 

The local surface heat flux, 7", was corrected for ra- 
diative heat transfer loss: 

q" = (k/ 1)(Tcu - T lc ) - «r(7£ - 2*,,) 

With c = 0.98, radiative losses were at most 6% of the 
net heat flux and typically much less. Variations in the 
radiative sink temperatures from T aw to between the exit 
adiabatic wall temperature and the inlet total tempera- 
ture produced no significant variations in Stanton num- 
bers. The Stanton number, St, was defined as follows: 

St = q" / pin Uin C r ( Tic - Taw) 

The local adiabatic wall temperature, T aw , is: 

Taw _ 1 ~ r 

r in “ 1 + 0.5(7- DA//N 

These definitions ensure that the measured Stanton num- 
ber is independent of liquid crystal temperature. The in- 
let total t emperature, T/ n , was measured by t wo probes 
located oik 1 axial chord upstream of the blades. The isen- 
tropic Mach number, A//*, was determined from CFD 
calculations, and a recovery factor of r — Pr 1/3 was used 
everywhere. 
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Fig. 4 Comparison of different blending approaches 
for turbulent eddy viscosity 

Uncertainty Analysis An uncertainty analysis was 
performed using the method of Kline and McClintock 
(11)53). The major sources of uncertainty are the con- 
ductance of the low conductivity layer, S(k/t)/(k/t) = 
±5%; the adiabatic wall temperature, ST a w = ±0 .()°C 
(±1.0°F); the liquid crystal temperature, STi c = ±0.3° C 
(db0.5°F); and the copper substrate temperature, STc u = 
±().^°( ' (±0.5° F). The uncertainty in T aw assumes no un- 
certainty in the recovery factor, r, although it is not pre- 
cisely known as shown by Schlichting (1979). The overall 
uncertainty in St was determined to be less than 13%. 
This maximum uncertainty is in regions of Stanton num- 
ber less than l x 10” 3 . For Stanton numbers greater 
than 2 x 10 -3 the uncertainty is less than 8%. Another 
possible source of bias uncertainty with this method is 
two-dimensional conduction within the low conductivity 
layer in regions of high spatial Stanton number gradients. 
A conduction analysis was performed to examine this er- 
ror. It was found that the distance from the minimum to 
t he maximum temperature in the region of largest spatial 
gradient was over 10 times the thickness of the low con- 
ductivity layer. A gradient of this magnitude was found 
to have virtually no adverse impact on the measurements. 

CFD ANALYSIS 

To illustrate where analysis improvements are 
needed, baseline CFD results were done using the 
three-dimensional Navier-Stokes analysis code, R\ C 3D, 
described by Chima and 5 okota (1990) and by 
Chima (1991). ('--type grids were generated using the 
method of Arnone ct al. (1992). Predictions were made 
using a 289 x 49 x 49 grid with a near wall spacing, yf < 1 . 
A uniform wall temperature was specified. 

A t wo- layer algebraic turbulence model, described by 
Cliinia et al. (1993). was used. The tavh blending be- 
tween the inner and outer layers was modified to give 
better agreement with the data. The turbulent eddy 
viscosity is given by: fit = F firo tan! \(^ri / ( 1 to)- Sub- 
scripts 77 and to refer to the inner and outer layers. In 
the unmodified blending F = 1. Since the tanh( 1) = 
0.70, this fi t is less than either the inner or outer value 
when both are equal. Other formulations, such as the 
Baldwin- Lomax model rely on an abrupt change between 
the in ner and outer regions. In the modifi ed blending 
F = >/(!.() + inin({fiT} / firo)’F (ftro / /'T/P)* Figure 4 



pitch 

Fig. 5(a) Inlet Mach numbers without turbulence grid 



pitch 


Fig. 5(b) inlet Mach numbers with turbulence grid 



Fig. 5(c) Inlet turbulence intensity with grid 



Fig. 5(d) Inlet turbulence length scale with grid 


compares the original and modified tank blending with 
an abrupt crossover model. The difference between the 
modified blending and the abrupt crossover model is small 
when firi < l { TO * the near wall region. However, the 
modified tank blending is higher after the crossover point, 
but rapidly returns to the outer layer value. 

Transition start was specified using Mayle's (1991) 
model. Tu at transition start was the upstream Tu times 
the ratio of the upstream velocity to the local isentropic 
velocity. The transition length model of Boyle and Si- 
mon (1998b), an extension of the one by Solomon et 
al. (1995) to include Mach number effects, was used. 

MEASUREMENT RESULTS 

Inlet Flow Measurements Inlet mean flow and tur- 
bulence quantities measured one axial chord upst ream of 
the blades are shown in Fig. 5. With no turbulence grid 
in place, the turbulence intensity was spatially uniform 
with values of 0.5% at /?c in = 0.5 x 10 (> and 0.25%' at 
Re in = 1.0 x 10 6 . 


NASA/TM—1 999-209296 


6 






z / span z / span 


M =0.98 

ex 



M = 1.32 

ex 



Fig. 6 Calculated isentropic surface Mach number 

Blade Static Pressures Figure 6 shows calculated 
isentropic surface Mach number distributions on the un- 
wrapped blade, since blade heat transfer is largely de- 
termined by flow features evident in the Mach num- 
ber distributions. The origin, s = 0, is the geomet- 
ric stagnation line at the design inlet flow angle. Giel 
et al., (1996a) showed that the Mach numbers obtained 
from the 3D Navier-Stokes code described earlier agreed 
well with blade loading measurements. The calculated 
values were used for heat transfer data reduction because 
they provided significantly better spatial resolution than 
was available from the experimental data. The surface 
Mach number distributions were not affected by Reynolds 
number variations, and their significance will be discussed 
in conjunction with the heat transfer results. 

Heat Transfer Measurements Stanton number con- 
tours on the unwrapped blade surface are shown in 
Figs. 7(a) through 7 ( /; ) * The eight flow conditions are 
given in Table 2. Exit Mach number variation effects are 
shown in alternate figures. Grid turbulence effects are 
seen by comparing alternate pairs of figures. Reynolds 
number effects are seen by comparing the first four fig- 
ures to the last four. Regions with no data available are 
shown as gray areas, and their extent depends on shadows 
and local heat transfer gradients. 

Figure 7(a) shows the Stanton number distribution at 
Rein — 0-5 x 10'\ M (T — 1.32, with the turbulence grid. 
The pressure side shows a relatively flat St distribution 
with a minimum above 3 x 10“ 3 , increasing to 5 x 10 -3 
near the trailing edge. The stagnation region is evident 
with a maximum St over 6 x 10" 3 . The Stanton number 
levels on the suction side of the leading edge initially drop 
to between 4 x 10“ 3 and 5 x 10 -3 , while those on the 
pressure surface decrease below 4 x 10” 3 . The grid causes 


Table 2 Description of blade heat 
transfer cases 


( ’ase 

Rf-in 

A/ej’ 

Tu grid 

1 

0.502 ±0.003 x 10 ,; 

1.322 ±0.003 

yes 

2 

0.500 ±0.002 x 10° 

0.985 ±0.001 

yes 

3 

0.498 ± 0.005 x 10 6 

1.313 ±0.003 

no 

4 

0.488 ±0.001 x 10 r ’ 

0.992 ±0.001 

no 

5 

1.003 ± 0.007 x 10 6 

1.322 ±0.001 

yes 

6 

1.000 ± 0.002 x 10 6 

0.981 ±0.001 

yes 

7 

0.999 ± 0.003 x 10 c 

1.314 ±0.001 

no 

8 

0.989 ± 0.003 x 10 r> 

0.998 ±0.001 

no 


(all repeatabilities based on 95% confidence limits) 


early transition, resulting in levels near midspan close to 
St = 5 x 10 -3 . The decelerating flow region (adverse 
pressure gradient) evident in Fig. 6 near s/s]>an — 0.5 
causes thickening of the turbulent boundary layer, result- 
ing in a minimum St of 3 x 10 -3 . From this point on. 
the effects of the secondary flows are evident . Many of 
these effects can be attributed to the vortex structures 
described by Goldstein and Spores ( 1 9 8 8 ) . The passage 
vortex and the pressure-side leg of the horseshoe vor- 
tex now approach the suction surface. Heat transfer is 
enhanced by the relatively cool secondary flow fluid im- 
pinging on the endwall regions of the suct ion surface. The 
vort ices lift off the endwall and approach midspan. Pre- 
vious measurements showed the vortices exit the rotor 
at approximately 25% span. The throat is located at 
s/span — 1.046 on the suction surface. The supersonic 
exit flow causes an oblique shock from the trailing edge 
of the adjacent blade to impinge on the suction surface. 
The effects of the shock/boundary layer interaction are 
evident at s/span = 1.3. The shock impingement causes 
a thickening of the boundary layer resulting in lower heat 
transfer. After the shock, secondary flows again cause 
higher heat transfer rates near the endwall. I he Stanton 
number then decreases as the suction surface boundary 
layer again grows as the flow approaches the trailing edge. 

Figure 7(b) shows the results for Case 2, which differs 
from the previous one by a reduction in exit Mach number 
from 1.32 to 0.98. The measurements show, as expected, 
that the heat transfer on the pressure surface and on the 
suction surface upstream of the throat are the same for 
these two cases to within the experimental uncertainty. 
Downstream of the throat, however, no shock is present 
to interrupt the suction surface boundary layer and the 
Stanton number contours continue smoothly to the trail- 
ing edge. The secondary flow effects on heat transfer are 
similar to those observed by Graziani et al. (1980) and by 
Blair (1994) in the region between the hub and midspan. 

Figure 7(c) shows measurements also taken at Rc m — 
0.5 x 10° and at M ex — 1.32, but with no inlet turbulence 
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Fig. 7 Measured Stanton number x 1 000 


grid. The absence of an upstream grid results in thicker 
endwall boundary layers, leading to more pronounced 
secondary flow effects near the pressure surface endwall. 
Away from the end walls, Stanton numbers are reduced by 
one-third to one-half. The pressure surface heat trans- 
fer is laminar-like, consistent with low turbulence level 
and strong favorable pressure gradients. As reported by 
Baughn et al. (1995), the pressure surface liquid crystal 
traces were far more nonuniform than those obtained with 
the grid. The traces had a finger-like appearance. 

The widt h and spanwise extent of the stagnation region 
contour line of St = (5 x 10“ 3 are both reduced over those 
of Case 1. The width reduction indicates a reduction 
in the peak heat transfer, as is expected with reduced 
freest ream turbulence. The reduction in spanwise extent 


results from the thicker inlet endwall boundary layers. 

No transition is observed between the stagnation re- 
gion and the adverse pressure region on the suction sur- 
face. The adverse pressure gradient causes tin' laminar 
boundary layer to grow rapidly and to possibly separate. 
The heat transfer reaches a minimum. St < 1 x J()~‘\ over 
a significant portion of the span. The induced freest roam 
turbulence from the secondary flows is apparently suffi- 
cient to trip the flow near the end walls thus preventing 
separation. The flow appears to separate, or to have the 
characteristics of a very abrupt transition. The transi- 
tion, or separation combined with a flow reattachment, 
causes an extremely high gradient in Stanton number. 
The levels increase from less than 1 x 10“ J to over 7 x 10“ 3 
in a stream wise distance less than one-tenth of the span. 
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Rivir et al. (1994) observed a similar sharp increase in 
heat transfer in a reattachment region downstream of a 
turbulent separation bubble. Further downstream, the 
Stanton number distribution of this no-grid case looks re- 
markably like that of the turbulence grid case, Case I 
Downstream of the shock, the levels are again generally 
in agreement with Case 1, but the distribution is much 
more distorted. Although the midspan symmetry appears 
to be quite good, the figure shows that two-dimensional 
midspan calculations would be inappropriate without ac- 
counting for midspan inflow. 

Figure 7(d) shows the measurement results for a sonic 
exit case. Case 4. Again, Re in — 0.5 x 10° with no grid 
present. Again, there is excellent agreement with the pre- 
vious supersonic exit case on the pressure surface and 
on the suction surface upstream of the throat. Down- 
stream of the throat, the contour lines generally follow 
the streamwise direction but again were more distorted 
than those of the turbulence grid case. 

Figure 7(c) shows the results of tin' first R(i t , = 
141 x 10 () case. Case 5. The turbulence grid was installed 
and the exit was supersonic. Comparing this case to 
Case 1 shows the effects of increased Reynolds number 
on heat transfer. The effects are small on the pressure 
surface because high freest ream turbulence gives early 
transition at either Reynolds number. As expected, the 
peak stagnation Stanton number is reduced about 20% 
with increased Reynolds number. Around midspan the 
high freestream turbulence and high Reynolds number 
cause very early suction surface transition, and St lev- 
els are hardly reduced from their peak stagnation region 
levels. The St levels decrease slightly as the suction sur- 
face boundary layer thickens. Because the flow is highly 
turbulent by this point, the slight adverse pressure gra- 
dient near s /span — 0.5 has almost no effect, on the heat 
transfer, and St levels remain at about 4 x 10 -3 until the 
shock impingement location. The core region following 
the shock impingement location again appears to be in- 
dependent of Reynolds number, t aking on values of about 
5 x lO" 3 . Near the endwall the Stanton number levels 
are nearly independent, of Reynolds number. Recall how- 
ever, that for equivalent Stanton numbers, doubling the 
Reynolds number doubles the heat transfer coefficient. 

The measurement results for this case clearly ver- 
ify that the blade design objective of more uniform heat 
transfer rates over the blade surface has been success- 
fully met . Those design objectives are most relevant for 
this case because the high Reynolds number and the high 
freestream turbulence levels are most similar to a high 
pressure turbine rotor operating environment. 

Comparing Fig. 7(f) to Fig. 7 ( r ) shows the effects of 
exit Mach number on the heat, transfer for two cases with 
Re: i„ = 14) x 1() G and with an inlet turbulence grid. 


pressure side suction side 



Fig. 8 Sample nondimensional surface temperature, 0 
(Re in = 0.5 x 1 0 6 ; M ex = 1 .32; no grid) 

The peak stagnation values for this case were just under 
5 x 1() -3 while for Case 5 they were just over 5 x 10 -3 . 
The two eases again look very similar on the pressure 
surface and upstream of the throat on the suction sur- 
face. Downstream of the throat, the contour lines again 
generally follow the streamwise direction. 

f igure '7(jt/) shows the results at R( j ti = 1.0 x 10*4 
M t r = 1.42, with no turbulence grid. The peak stag- 
nation region Stanton number is just under 4 x 10 -3 . 
The pressure surface is transitional because of the high 
Reynolds number. Suction surface transition is more 
complicated than any of the previous cases. At about 
40%, 50%, and 70% span, the flow starts to transition 
just, downstream of the stagnation region and continues 
the transition process relatively slowly. At 40%. and ()0% 
span however, the flow remains laminar until it reaches 
the slight adverse pressure gradient at s/span = 0.54 
where it appears to separate, then reattach and become 
fully turbulent. The behavior at these two discrete span- 
wise locations resembles that of the entire core region of 
t lie R( in — 0.5 x 10*4 no grid cases. The most likely 
cause of the relatively early transition spots was the dots 
painted on the surface at 44%., 50%, and 07% span at 
s/spau = 0.44. The dots were barely perceptible to the 
touch, but the high Reynolds numbers apparently made 
the boundary layer very sensitive to small disturbances. 
Following this region, the Stanton number contours again 
resemble all of the other supersonic exit cases. 

Figure 7(h) shows that the pressure surface and the 
suction surface upstream of the throat again agree very 
well with the previous supersonic* exit case. Downstream 
of the throat, the Stanton number contours again follow 
the streamwise direction. 

The measurement technique produces a blade sur- 
face that is not uniform in either surface temperature or 
surface heat flux. A sample non-dimensional blade sur- 
face temperature, 0 , is shown in Fig. 8 for the Case 4 
which had the largest 0 variation. 0 is independent of 
the particular measurement, temperatures and is inversely 
proportional to the Stanton number, so that boundary 
conditions for all cases can be found from the Stanton 
numbers. 
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(d) Re jn = 0.5 X 1 0 6 ; M ex = 0.98; no grid (h) Re* = 1 .0 x 1 0 6 ; M ex = 0.98; no grid 


Fig. 9 Measured and predicted Stanton no. x 1000 

An example using Case 3 (Fig. 7c) will now he dis- approximately T sur j act /T{ n — 1.73, showing that t he 
cussed to relate actual temperatures to values of 0. Mea- present measurement technique produces a thermal 
surements near the stagnation region with Ti c ^/T' n — boundary condition much closer to a uniform surface tem- 
1.05) showed St =6.1 x 1U“ 3 and 0 — 1.09. These mea- perature than a uniform heat flux. 

surements required a copper temperature of T Cu /T' n = The liquid crystal Stanton number data were int.er- 

1.18. Because t he copper temperature was essentially uni- polated onto lines at z/span — 0.10, 0.25, and 0.50. Line 
form, regions of low heat transfer had higher surface tern- plots of this data are shown in Figs. 9(a) through 9(fc). 
peratures. For example, in the separated flow region the The contour plots of Fig. 7 show overall trends and span- 
minimum Stanton number was approximately 0.7 x 10 -a . wise variations, but line plots can be useful for extracting 
6 here was approximately 8.6, resulting in a surface tern- more details of the data. Computational results at the 
perature of T sur f ace /Tj n — 1.16. Had these measure- same spanwise locations are also shown in the figures, 
ments been made with a uniform heat flux technique, the 
surface temperature in the low St region would have been 
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(a) Re in = 0.5 x 10 6 ; M ej< = 1.32; with grid (e) Re in = 1.0 x 10 6 ; M ex = 1.32; with grid 





(d) Re in = 0.5 x 1 0 6 ; M ex = 0.98; no grid (h) Re* = 1 .0 x 1 0 6 ; M ex = 0.98; no grid 

Fig. 1 0 Predicted Stanton number x 1000 


Some general comments rail be made that are com- 
mon to all of the cases presented above. First, the peak 
stagnat ion heat transfer occurs to the suction side of the 
geometric stagnation line, in contrast to the aerodynamic 
stagnation line which is slightly towards the pressure sur- 
face. Second. Stanton numbers of 3 x 1()“ 3 were observed 
in the shock impingement region of all supersonic exit 
cases, regardless of Reynolds number or turbulence grid. 
Third, secondary flows significantly increase suction sur- 
face heat transfer rates near the endwalls. Fourth, the 
data with the grid in place was smoother than the no-grid 
data, which was felt to be due to the strong mixing effects 
of the freest ream turbulence. Finally, excellent midspan 
symmetry is observed. 


COMPUTATIONAL HEAT TRANSFER 

The predicted heat transfer rates are shown in 
Figs. 10(a) through 10(A) in the same order as the exper- 
imental data. The analysis assumed midspan symmetry, 
but full span solutions are shown for comparisons with 
data. An inlet Tu — 8% was used for predictions with 
the grid installed, and Tu = 1% for the no grid cases. 

Low Reynolds number comparisons Figure 10(a) 
shows results for Rc in - 0.5 x 10 6 , M ex = 1.32, and with 
the grid. In the analysis, the turbulence intensity affected 
start of transition and gave increased laminar region heat 
transfer. Leading edge and nearby suction surface heat 
transfer are well predicted. The analysis and data both 
show suction surface minimums near s/span — 0.5. The 
minimum experimental Stanton number in this region is 
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3 x 10” 3 . At ti/span ~ 0.5 predictions show a region at 
this level, and a minimum at 2.5 x 10“ 3 close to one quar- 
ter span. The midspan peak Stanton number of 4 x 10 -3 
is downstream of the experimental peak, indicating the 
transition length is over predicted. However, overall the 
agreement with the data in the transition region is good. 
At s/span = 0.5 a high heat, transfer region forms near 
the end wall. This region expands downstream consistent 
with the data. The peak experimental Stanton number 
is 7 x 10 -3 , while the peak predicted value is 0.5 x 10” 3 . 
Heat transfer behavior in t his region is due to secondary 
flows as discussed by Goldstein et al. (1995), and is not 
just due to (toss channel flows passing over the unheated 
endwall. An analysis with the blade and endwall both 
heated showed blade heat transfer very similar to the pre- 
dictions. The effect of the shock on the heat transfer is 
clearly seen, as it was in the data. After the shock, the 
predicted midspan heat t ransfer decreases, but this is not 
seen in the data. Close to the end walls both the analysis 
and t he data show regions of increased heat transfer af- 
ter the shock. The peak computed value in this region is 
1.5 x nr 3 , while the peak data value was 5 x 1U“ 3 . 

Both the computations and experiment show a very 
rapid decrease in pressure surface heat t ransfer away from 
the leading edge. The minimum experimental contour 
level is 4 x 1 ()~ 3 , and encompasses nearly two thirds of 
the pressure surface. The absolute minimum measured 
pressure surface Stanton number was close to A x H) -3 . 
Except, close to the endwall, the minimum predicted Stan- 
ton number is 2.5 x 10 -3 . Transition was predicted to 
begin at close to A0% of the surface distance, and was 
not complet e until after 50% of the surface distance. Un- 
derpredicting the minimum Stanton number appears to 
imply a poor start of transition prediction. However, if 
the analysis underpredicted the effects of freestream tur- 
bulence on the pressure surface heat transfer, it. would 
also under predict- the heat transfer, even if the start and 
length of transition were correctly predicted. Close to the 
endwall at 40% of the pressure surface distance a region of 
low heat transfer is predicted, but is not seen in the data. 
At (dose to 90% of the surface distance there are Stanton 
number levels of 5 x 10~ 3 in both the experimental and 
computational contours. 

Figure 10(6) shows predictions for the lower exit 
Mach number case. A change in heat transfer is only 
seen from the t hroat region to the trailing edge of the 
suction surface. At this lower Mach number there is lit- 
tle evidence that the much weaker shock affects the heat, 
transfer. Over the last 4i)% of the suction surface the 
Stanton number lines are nearly horizontal. This loca- 
tion is governed by the cascade aspect ratio. The cas- 
cade symmetry forces the secondary and vortex flows to 
stop moving closer to midspan at some chordwise loca- 


tion. Close to midspan the analysis shows lower than 
measured heat transfer. In addition to a deficiency in the 
turbulence modeling, the analysis could be overpredicting 
the amount of secondary flows forced into the midspan 
region, giving a thicker boundary layer and reduced heat 
transfer. 

Figure 10(c) shows heat transfer levels on the pres- 
sure surface, and the suction surface prior to transition, 
are affected when the grid is absent. The pressure surface 
heat transfer is much lower at the 1%; inlet Tn level. At 
the leading edge the midspan Stanton number is slightly 
less than 4.5 x 10“ 3 . This corresponds to a Frossling num- 
ber, Nu d /y/Rc d < of 1.04. The experimental leading edge 
Frossling number is higher, 1.39. There is no mechanism 
for the model to augment the heat transfer at low tur- 
bulence intensity, so the analysis gives the same result as 
laminar flow at the stagnation point of a cylinder. The 
calculations predicted laminar flow on t he pressure sur- 
face. In the experimental data there is a large pressure 
surface area enclosed within the St = 2x 10~ 3 contour. In 
the calculations there is a minimum contour of 1.5 x 10~ 3 
and at midspan it extends into the region where the ex- 
perimental data show increased heat transfer, ('loser to 
the pressure surface trailing edge the experimental Stan- 
ton number level is 3 x 10 -3 , while the calculated value 
is only 2 x 10 -3 . These results indicate that there is a 
mechanism increasing the pressure surface heat transfer 
that the analysis does not account for. 

Suction surface heat transfer also differs bet ween the 
high and low Ta cases due to turbulence modifying t he 
transition start location. Data show a minimum Stan- 
ton number contour of J x 10~ 3 . The prediction shows a 
minimum contour level of 1.5 x 10“ 3 . enclosing a some- 
what larger area. The start of transition is accurately 
predicted. Prior to transition, the suction surface of the 
blade show’s the same relative behavior as the pressure 
surface. In these laminar regions the predicted heat, trans- 
fer is low'er than the data. 

In the midspan region, after the minimum suction 
surface heat transfer, the predicted heat transfer shows a 
much lower peak than the data. These differences may 
be due to flow' separation and reattachment , rather than 
transition. The analysis predicted a low' shear value at 
the start of transition, and a very short transition region 
caused by the adverse pressure gradient. If the actual 
transition w r as slightly later than predicted, separation 
would have occurred. Very high heat transfer is expected 
at reattachment. The experimental heat transfer show's 
two small contour areas above and below midspan at. a 
level of St = 7 x 10 -3 . The analysis showed a con tom- 
level in the same region between 3.5 x 10“ 3 and 4.0 x 10~ 3 . 
The algebraic turbulence model used did not, correctly 
predict the heat transfer in this region, and improved 
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turbulence models should give better agreement. Cal- 
culations were done with the wall temperatures varying 
in the same manner as that shown in Fig. 8, but did not 
improve the agreement with data. 

The increased near-endwall heat transfer due to sec- 
ondary flows is very similar to the no grid case, and is 
well predicted. The peak predicted Stanton number is 
6.5 x 1 0 — 3 , while the peak measured value is 7 x 10“ 3 . 
The behavior in the shock impingnient region is also well 
predicted. Both data and predictions show contour levels 
of x 10 -3 . Downstream of the shock, near midspan, the 
data show Stanton number levels above 4 x 1()“ 3 . The 
prediction shows a decrease in heat transfer towards the 
trailing edge. This behavior is similar to that for the cor- 
responding high Tv case. Away from midspan the grid 
case data show a higher average heat transfer level, and a 
more chaotic distribution. The analysis, however, shows 
similar behavior but at a slightly lower level. The surface 
pressure distributions were the same for the high and low 
Tit analyses. There was no mechanism in the analysis 
that would account for higher freest ream turbulence af- 
ter the flow was fully turbulent . 

Figure 10(d) shows that decreasing Mach number 
from 1.32 to 0.98 for the no grid case resulted in changes 
similar to those for the same Mach number change at the 
high Tu grid case. 

High Reynolds number comparisons The primary 
effect of increasing the Reynolds number by a factor of 
two on the heat transfer distribution is to change the tran- 
sition location. Prior to transition the Stanton number is 
expected to decrease by \/2 when the Reynolds number is 
doubled. For fully turbulent flow the Stanton number is 
expected to decrease by 13% when the Reynolds number 
is doubled, assuming St (x Re~ [) -. 

The first comparison is between the two Reynolds 
numbers at A/ ex = 1.32 and the grid installed. Fig- 
ure 10(c) shows that downstream of the leading edge on 
the suction surface the Stanton number is higher for the 
higher Reynolds number, as it is in the data. This is 
due to transition moving forward at the higher Re. The 
leading edge heat transfer decreases only by about 15%. 
However, this is consistent with the high turbulence in- 
tensity and high Red . The correlation of Van Fossen 
et al. (1995) predicts a 2(5% reduction in St due to the 
change in Reynolds number. However, the correlation 
also predicts that at higher Frossling numbers, seen in 
the data, the effect of a Reynolds number variation is 
less. Elsewhere on the blade, the change in Stanton num- 
ber is as expected. The flow is turbulent, and there is 
a decrease of 10-20%. in the Stanton number when the 
Reynolds number is doubled. 


Comparing Figs. 10(c) and 10(/) shows the same heat 
transfer changes at the high Reynolds number for a Mach 
number change as were seen at the low Reynolds number. 
Overall, the agreement with the data is good. However, 
again the predictions show too low a heat transfer on the 
suction surface midspan region downstream of the throat. 

Figure 10(g) shows that, with no grid in place, in- 
stead of decreased St when Re doubles, the St anton num- 
ber more than doubles near the pressure surface trailing 
edge. This is due to transition. At the lower Reynolds 
number, pressure surface transition did not, occur. At the 
higher Reynolds number transition was predicted to be- 
gin near 30% of pressure surface distance, and end at 60% 
of the surface distance. This is consistent with t he experi- 
mental dat a, although the change in St is not as dramatic 
in the data. Suction surface transition was predicted to 
occur at nearly the same location for both Reynolds num- 
bers. Prior to transition the lower Reynolds number had 
a region between St = 1.5 and 2 x 10" 3 . For the higher 
Reynolds number the same region shows Stanton num- 
bers less than 1.5 x 10 -3 . This is consistent with what 
one would expect for laminar flow. Figure 10 (It) shows 
that the change in Stanton number with decreased Mach 
number is the same at both Reynolds numbers. 

Overall the analysis underpredicts the heat transfer, 
which may be partly due to differences in the manner 
in which the heat transfer coefficient was defined. For 
consistency with the experimental data the Stanton num- 
bers are based on the wall heat flux divided by the dif- 
ference between the wall temperature and the adiabatic 
wall temperature. The same adiabatic wall temperature 
distribution was used to determine the experimental and 
computational Stanton numbers. Therefore, differences 
are due to differences in the wall heat flux. Calculations 
were done with different wall temperatures. From the 
resultant heat fluxes the adiabatic wall temperature was 
calculated for each grid point on the blade surface. The 
recovery factor determined from the Navier-Stokes anal- 
ysis was closer to unity than the recovery factor used to 
calculate the Stanton numbers. The recovery factor was 
also greater than most of the recovery factors given by 
Schlichting (1979). If the inlet total temperature was 
used to calculate Stanton numbers instead of the local 
adiabatic wall temperature, predicted Stanton numbers 
would be up to 40% higher. Differences arising from dif- 
ferences in the recovery factor only help to explain dif- 
ferences in the suction surface heat transfer. Only heat 
transfer in high Mach number regions is affected by vari- 
ations in the recovery factor. 
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SUMMARY AND CONCLUSIONS 

The primary effect of a factor of two variation in 
Reynolds number was to move the transition locations. 
For the low inlet turbulence intensity cases, flow on the 
suction surface remained laminar up to an adverse pres- 
sure gradient region. Subsequently, the Stanton number 
rose rapidly to a high level. At low Tu, pressure surface 
transition was observed at high, but not at low, Reynolds 
numbers. The shock in the supersonic exit cases impinged 
on the suction surface, where measurements and calcu- 
lations showed St anton number levels of 3 x 10 -3 inde- 
pendent of Reynolds number and turbulence grid. The 
subsonic cases showed Stanton number contours gener- 
ally following the expected streamlines and remaining 
unchanged after the throat. The data showed that the 
turbulence grid increased leading edge heat transfer and 
moved the transition locations forward. For the high 
Reynolds number and turbulence grid case the suction 
surface transition moved almost to the stagnation point. 

Kven with high turbulence the midspan leading edge 
heat transfer did not greatly exceed that seen elsewhere 
along midspan. The good spatial resolution due to the 
large scale and the liquid crystal technique allowed the 
secondary How effects to be clearly quantified. 

Comparing computational and experimental results 
illust rated regions of good agreement and regions where 
modeling improvements are needed. Transition was well 
predicted, as were the effects of secondary flows on the 
suction surface heat transfer. High heat transfer in the 
near end wall region was well predicted. The leading edge 
heat transfer was well predicted for the high turbulence 
intensity rases, but was underpredicted for the low tur- 
bulence intensity cases. For the low Reynolds number, no 
grid rases, the analysis predicted a very short suction sur- 
face transition region, due to an adverse pressure gradi- 
ent. Experimentally, the laminar flow may have separated 
and reat tached rather than smoothly transitioned, which 
may be the reason that the analysis underpredicted the 
heat transfer. The analysis also underpredicted the heat 
t ransfer near midspan downstream of the throat on the 
suction surface. This could be due to turbulence model 
deficiencies, or an overprediction of the amount of fluid 
forced into the midspan region by secondary flows. The 
latter would cause the calculated midspan boundary layer 
to be thicker than the experimental one, leading to lower 
heat transfer. 

The blade heat transfer data presented here, the end- 
wall heat transfer data of (tie I et al. (1998), along with 
the aerodynamic data of Ciel et al. (1996a) comprise a 
complete set of data for CFD code and model validation. 


Electronic data tabulations for all eight cases, including 
raw digitized data to allow contour plotting at any level, 
are available upon request. 


REFERENCES 

Arnone, A, Liou, M.-S. and Povinelli, L. A., 1992, 
“Navier- Stokes Solution of Transonic Cascade Flows Using 
Non-Periodic C-Type Grids," AIAA Journal of Propulsion 
and Poiver. Yol. 8, No. 2, pp. 410-417. 

Arts, T., Duboue, J.-M., and Rollin, G., 1997, “Aero- 
Thermal Performance Measurements and Analysis of a Two- 
Dimensional High Turning Rotor Blade," ASME paper 97- 
GT- 120. 

Baughn, J. W., Butler. R. .J.. Byerley, A. R., and Rivir, 
R. B., 1995, u An Experimental Investigation of Heat Transfer, 
Transition and Separation on Turbine Blades at Low Reynolds 
Number and High Turbulence Intensity," ASME paper 95- 
YVA/HT-25. 

Blair, M. F.. 1994, “An Experimental Study of Heat 
Transfer in a Large-Scale Turbine Rotor Passage," ASME 
Journal of Turbomachinery, Vol. 116, No. 1. pp. 1-13. 

Boyle, R. J., Lucci, B. L.. VerhofF, V. G., Camperchioli, 
YV. P. . and La, H., 1998a, “Aerodynamics of a Transitioning 
Turbine Stator Over a Range of Reynolds Numbers," ASME 
paper 98-GT-285. 

Boyle, R. J,, and Simon, F. F.. 1998b, “Mach Number Ef- 
fects on Turbine Blade Transition Length Prediction," ASME 
paper 98-GT-367. 

Chima, R. V., and Yokota, J. YV., 1990, “Numerical 
Analysis of Three-Dimensional Viscous Internal Flows," AIAA 
Journal, Vol. 28, No. 5, pp. 798-806. 

Chima, R. V,, 1991, “Viscous Three- Dimensional Calcu- 
lations of Transonic Fan Performance," AGARD Propulsion 
and Energetics Symposium on Computational Fluid Mechan- 
ics for Propulsion, San Antonio, Texas, May 27-31. 

Chima, R. V., Giel, P. YV., and Bovle, R. J., 1993. “An 
Algebraic Turbulence Model for Three-Dimensional Viscous 
Flows," AAIA paper 93-0083, (NASA TM- 105931). 

Civinskas, K. C., Boyle. R. J., and McConnaughey, H. V., 
1990, “Turbine Blading Designed for High Heat Load Space 
Propulsion Applications." AIAA Journal of Propulsion and 
Power, Vol. 6, pp. 598-611. 

Dunn, M. G., Kim, J., Civinskas, K. C., and Boyle, R. J., 
1994, “Time- Averaged Heat Transfer and Pressure Measure- 
ments and Comparison with Predictions for a Two-Stage Tur- 
bine," ASME Journal of Turbomachinery . Vol. 116, pp. 14-23. 


N AS A/TM— 1 999-209296 


14 



Giel, P. W., Thurman, D. R., Lopez, I., Bovle, R. J., 
Van Fossen, CL J., Jett, T. J., Camperchioli, W. P., and 
La. H.. 1996a, ‘‘Three-Dimensional Flow Field Measurements 
in a Transonic Turbine Cascade," ASME paper 96-GT-113. 
presented at the ASME International Gas Turbine Confer- 
ence, Birmingham, England, June 10-13. 

Giel, P. W.. Sirbaugh, J. R., Lopez, L, and Van Fos- 
sen. G. J.. 1996b, “Three-Dimensional Navier-Stokes Analysis 
and Redesign of an Imbedded Bellmouth Nozzle in a Turbine 
Cascade Inlet Section." ASME Journal of Turbornachinery , 
Vol. 118, No. 3, pp. 529-535. 

Giel, P. W., Thurman, D. R., Van Fossen. G. J., Hippen- 
steele, S. A., and Bovle, R. J., 1998, “Endwall Heat Transfer 
Measurements in a Transonic Turbine Cascade," ASME Jour- 
nal of Turbornachinery. Vol. 120, No. 2, pp. 305-313. 

Goldstein, R. J.. and Spores. R. A., 1988. “Turbulent 
Transport on the Endwall in the Region Between Adjacent 
Turbine Blades," ASME Journal of Heat Transfer , Vol. 110. 
No. 4. pp. 862-869. 

Goldstein. R. J., Wang, H. P.. and .Jabbari, M. V.. 1995, 
“ The Influence of Secondary Flows Near the Endwall and 
Boundary Layer Disturbance on Convective 'Transport From a 
Turbine Blade," ASME Journal of I urbornachin e ■ ry. \ol. 11/ , 
No. 4, pp. 657-665. 

Graziani. R. A.. Blair, M. F., Taylor. R. J.. and Mayle, 
R. E., 1980, “An Experimental Study of Endwall and Airfoil 
Surface Heat Transfer in a Large Scale Turbine Blade Cas- 
cade," ASME Journal of Engineering for Power. Vol. 102, 
No. 2, pp. 1-11. 

Guenette. G. R., Epstein. A. H., Giles, M. B., Haimes, R., 
and Norton, R. J. G.. 1989, “Fully Scaled Transonic 'Turbine 
Rotor Heat Transfer Measurements," ASME Journal of Tur- 
bomachinery. Vol. Ill, No. 1, pp. 1-7. 

Hippensteele, S. A., Russell. L. M., and Torres, F. J., 
1988, “Local Heat -'Transfer Measurements on a Large Scale- 
Model Turbine Blade Airfoil Using a Composite of a Heat er 
Element and Liquid Crystals," ASME! Journal of Engineering 
for Gas Turbines and Power. Vol. 107, pp. 953-960. 


Joslyn, D., and Dring, R., 1992, “Three-Dimensional 
Flow in an Axial Turbine: Part 1 - Aerodynamic Mecha- 

nisms," ASME Journal of Tutbomac hinery. Vol. 114, No. L 
pp. 01-70. 

Kline, S. J.. and McOlintock, F. A., 1953, “Describing 
Uncertainty in Single-Sample Experiments," Mechanical En- 
gineering, Vol. 75, Jan., pp. 3-8. 

Martinez-Bot as, R. F.. Lock, G. D., and Jones, T. V ., 

1994, “Heat Transfer Measurements in an Annular Cascade 
of Transonic Gas Turbine Blades Using the Transient Liquid 
Crystal Technique," ASME 94-GT-172. 

Mavle, R.E., 1991, “The Role of Laminar- Turbulent 'Tran- 
sition in Gas l urbine Engines," ASME Journal of Turboma- 
chinery. Vol. 113, pp. 509-537. 

Moffat, R. J., 1990. “Experimental Heat 'Transfer." Proc. 
of tin Ninth hit 'I Heat Transfer Conf., Jerusalem. Israel, 
Vol. 1, pp. 187-205. 

Rivir, R. B., Johnston, J. P., and Eat-on, J. K., 1994, 
“Heat Transfer on a Flat Surface Under a Region of Turbu- 
lent Separation." ASME] Journal of Turbomachinery. Vol. 1 16, 
No. 1, pp. 57-02. 

Schlicht ing, IF, 1979, Boundary- Layer Theory , Seventh 
Edition. McGraw-Hill, New York, p. 714. 

Solomon, YY. J.. Walker, G. J., and Gostelow, J. P. . 

1995. “Transition Length Prediction For Flows With Rapidly 
(■hanging Pressure Gradients," ASMF] paper 95-GT-24L 

Timlin. R. I).. Howe, D. <\, and Singer, I. I)., 1982. “En- 
ergy Efficient. Engine High-Pressure Turbine Detailed Design 
Report." NASA CR- 165608. 

Van Fossen. G. J.. Simoneau, R. L, and Oiling, (\ Y ., 
1995, “Influence of Turbulence Parameters, Reynolds Number, 
and Body Shape on Stagnation-Region Heat 'Transfer," ASME] 
Journal of Heat Transfer. Vol. ID, pp. 59 <-603. 

Verhoff, V. G.. Camperchioli, YV. P., and Lopez, L. 1992, 
“Transonic Turbine Blade Cascade Testing Facility." A1AA 
Paper No. 92-4034, NASA TM- 105646. 


15 


NASA/TM— 1 999-209296 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway Suite 1204. Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0168), Washington, DC 20503. 


1. AGENCY USE ONLY {Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 


August 1999 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


WU-523-26- 1 3— (K) 
1L 161 102AH45 
NAS3-98008 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


Blade Heat Transfer Measurements and Predictions in a Transonic 
Turbine Cascade 


6. AUTHOR(S) 

PW. Giel. G.J. Van Fossen, R.J. Boyle, D.R. Thurman and K.C. Civinskas 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Glenn Research Center 
Cleveland. Ohio 44135-3191 
and 

U.S. Army Research Laboratory 
Cleveland. Ohio 44135-3191 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 

Washington. DC 20546-(XX)l 

and 

11, S. Army Research Laboratory 
Adelphi. Maryland 20783-1 145 


11. SUPPLEMENTARY NOTES 

Prepared for the 1999 International Gas Turbine and Aeroengine Congress cosponsored by the American Society of 
Mechanical Engineers and the International Gas Turbine Institute, Indianapolis, Indiana, June 7-10, 1999. P.W. Giel, 
Dynacs Engineering Company, Inc., Brook Park, Ohio; (work funded under NASA Contract NAS3-98008); G.J. 
VanFossen and R.J. Boyle. NASA Glenn Research Center; D.R. Thurman and R.C. Civinskas, U.S. Army Research 
Laboratory, NASA Glenn Research Center. Responsible person, R.J. Boyle, organization code 5820, (216) 433-5889. 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 1 12b. DISTRIBUTION CODE 


E- 1 1775 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM— 1999-209296 
ARL-TR-2029 


Unclassified - Unlimited 
Subject Categories: 02, 07 and 34 


Distribution: Nonstandard 


This publication is available from the NASA Center for AeroSpace Information, (301 ) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

Detailed heat transfer measurements and predictions are given for a turbine rotor with 136° of turning and an axial chord 
of 12.7 cm. Data were obtained for inlet Reynolds numbers of 0.5 and l.OxlO 6 , for isentropic exit Mach numbers of 1.0 
and 1 .3, and for inlet turbulence intensities of 0.25% and 7.0%. Measurements were made in a linear cascade having a 
highly three-dimensional How field resulting from thick inlet boundary layers. The purpose of the work is to provide 
benchmark quality data for three-dimensional CFD code and model verification. Data were obtained by a steady-state 
technique using a heated, isothermal blade. Heat fluxes were determined from a calibrated resistance layer in conjunction 
with a surface temperature measured by calibrated liquid crystals. The results show the effects of strong secondary 
vortical flows, laminar-to-lurbulent transition, shock impingement, and increased inlet turbulence on the surface heat 
transfer. 


14. SUBJECT TERMS 

Gas turbines; Heat transfer; Transonic flow; Supersonic turbines; Cascade wind tunnel; 
Liquid crystals; Turbomachinery; Turbulent heat transfer; Convective heat transfer; 
Stanton number 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


15. NUMBER OF PAGES 
21 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 














